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Abstract 

We consider the possibility that sterile neutrinos exist and admit a description as unparticles; 
we call these unsterile neutrinos. We then examine the nature of neutrino oscillations in a model 
where an unsterile can mix with two active flavors with a very simple mass matrix of the see-saw 
type. Despite these simplifications, we find a number of remarkable features, all of which will 
occur generically when more realistic cases are considered. These include momentum dependent 
mixing angles, "invisible" decay widths for the unsterile-like mode, as well as the inheritance of a 
non- vanishing spectral density for the massive active-like modes. We also obtain the disappearance 
and appearance probabilities for the active-like neutrinos and find remarkable interference effects 
between the active and unsterile neutrinos as a consequence of threshold effects, yielding new 
oscillatory contributions with different oscillation lengths. These effects are only measurable on 
short baseline experiments because there both probabilities are suppressed as compared to mixing 
with a canonical sterile neutrino, thereby providing a diagnostics tool to discriminate unsterile 
from canonical sterile neutrinos. We conclude with a discussion of whether these new phenomena 
could aid in the reconciliation of the LSND and MiniBooNE results. 
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I. INTRODUCTION 



Neutrino masses and oscillations are now an established phenomenon and indisputable ev- 
idence of physics beyond the standard model. Although the origin and scale of these masses 
remains a challenging question, the see-saw mechanism provides a compelling framework on 
how to explain the small active neutrino masses [1]. 

One particular extension of the Standard Model consists of adding the so called sterile 
neutrinos, which are massive SU(2) singlets. Scenarios involving sterile neutrinos with 
mass in the GeV range have been proposed as explanations of two seemingly unrelated 
and unexpected phenomena: an excess of air shower events at the SHALON gamma ray 
telescope [2], and the prominent peak of electron- neutrino events above background for 
300 MeV < Ey < 475 MeV at the MiniBooNE [3] . Examples of such scenarios include a 
3+2 oscillation scheme involving three active neutrinos and two sterile neutrinos proposed 
in Ref. [1], and an alternative explanation invoking the radiative decay of the heavy sterile 
neutrino with a small magnetic moment proposed in Ref. [3]. 

Recently, Georgi [6] suggested another extension of the Standard Model, which includes a 
"hidden" conformal invariant sector with a non-trivial infrared fixed point. This conformal 
sector can be realized by a renormalization flow toward the infrared through dimensional 
transmutation [71 [8]. Below the dimensional transmutation scale, there emerges an effec- 
tive interpolating field, the unparticle field, that features an anomalous scaling dimension. 
Various studies recognized important phenomenological lUHTB] , astrophysical [T7HTU] and 
cosmological [20H28] consequences of unparticles. The consequences of the existence of scalar 
unparticles in the neutrino sector have also been studied. A scenario where the heaviest mass 
eigenstates corresponding to the mixed — z/g decaying into the lightest eigenstate and a 
scalar unparticle has been suggested as an explanation of the MiniBooNE anomaly [29l 130] . 
Furthermore, unparticle contributions to the neutrino-nucleon cross section and its influence 
on the neutrino flux expected in a neutrino telescope such as IceCube have been reported 
in Ref. [31]. 

The question arises whether interesting physics might be found by combining these two 
ideas. Thus, suppose that sterile neutrinos do exist, and furthermore, they admit a descrip- 
tion in terms of unparticles. What consequences arise for neutrino physics, and in particular, 
neutrino oscillations? 
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In Ref. [32], we considered the sterile neutrino to be an unparticle, namely an unsterile 
neutrino, and studied the mixing of an active and an unsterile neutrino, which revealed 
unexpected momentum dependence of the mixing angles. Here, we extend that model by 
including a second active neutrino, and study the consequences in neutrino oscillations. To 
be precise, we consider a 2 + 1 scheme in which two massless active neutrinos mix with 
an unsterile neutrino via a simple mass matrix of the see-saw type that mixes the active 
neutrinos only with the unsterile one. Despite the simplicity of the mass matrix, a number 
of results arise that will persist in more complex models. 

We find that the unsterile-active mixing angle depends on the four-momentum and that 
the propagating modes consist of a massless active-like mode, a massive active-like mode 
and an unsterile-like mode. The structure of propagators of the massive active-like mode 
and the unsterile-like mode are similar to the ones found in [32] . 

To study the dynamics of active-active oscillation, we generalize the conventional quan- 
tum mechanical description of neutrino oscillations in terms of the propagators of mass 
eigenstates and calculate the appearance and disappearance probabilities. The oscillation 
dynamics reveals remarkable interference phenomena associated with the threshold of the 
inherited spectral density for the massive active-like and the complex unsterile-like pole. 
The momentum dependent mixing angles along with these novel interference contributions 
yield new oscillatory contributions to the appearance and disappearance probabilities of the 
active neutrinos with a different oscillation length. These new contributions will be mani- 
fest as quantum beats as a function of energy and baseline in oscillation experiments. The 
non-oscillatory contributions to these probabilities are suppressed with respect to the case 
of a canonical sterile neutrino, which has no anomalous dimension. 

Recently, in Ref. |33], it has been proposed that energy dependent mixing and or oscil- 
lation parameters may reconcile the LSND [M] and MiniBooNE [3] data. This suggestion 
hinges on the possibility of "exotic" sterile neutrinos whose mixing with the active ones 
features novel energy dependent mixing angles and/or oscillation parameters. The energy 
dependence may reconcile the results of both experiments which are performed in different 
energy regimes (~ 40MeV for LSND, vs. ~ GeV for MiniBooNE). Given the energy depen- 
dence we find in our analysis, we were motivated to see whether our model could be of use 
in understanding how the LSND and MiniBooNE data mesh together. For our model, we 
find that we cannot reconcile the data, but that might be attributable to not having enough 
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structure in the mass matrix. 



II. 2 + i SCHEME 



The Lagrangian density for the unsterile in momentum space is 
where 



Fip) 



+ M2 - iO+ 
A2 



0<T] <1. 



(2.1) 



(2.2) 



A is the scale below which the low energy dynamics is dominated by the infrared fixed point 
of the conformal sector. Below this scale the unparticle is described by an interpolating 
field whose two point correlation function scales with an anomalous dimension. Consistency 
of the unparticle interpretation requires that M < A, where M is the infrared cut-off or 
unparticle threshold. This Lagrangian density can be understood using a renormalization 
group resummation argument discussed in Ref. [3^ and references therein. 

We consider a see-saw type mixing with two active Dirac neutrinos, such that the full 
momentum space Lagrangian density is given by 

/ u 



C 



VaA~P) VaA-p) i'ui-p) 







(2.3) 



—mi 

^ -m2 i^a2{p) 

y-mi -1712 {^-M)F{p) j y^c/(p) j 

with 2 &re the active neutrinos and we assume mi^2 ^ M < A leading to a see-saw 
hierarchy of masses. Following [32], let us introduce a rescaled unsterile field 

vu = VfIp) iJu, (2.4) 
such that the Lagrangian density is now given by 

C={u„^ Vu ) (^I - M) 



■'02 



(2.5) 



where I is the identity in the "flavor" space, and the "mass" matrix is 



M 







ibll \ 



^ 



(2.6) 
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We note that M has a zero eigenvalue, which aids considerably in the diagonalization. We 
reiterate that this mass matrix describes a simplified model, chosen only so as to exhibit 
the mixing that we want to describe (active-unsterile) while suppressing effects we consider 
extraneous to the discussion. 

This mass matrix can be diagonalized exactly and the fields corresponding to the propa- 
gating or "mass" eigenstates are given by 

-1712 J^ai + mi Va2 7^ 
Vi — , [^Z.i) 

m 



2mi Ua, + 2m2 z/^^ - (VM^F + Am"^ - M^fP 

^2 = ^ -d ' (2-8) 

\2M'^F + Sm? — 

2mi Ua, + 2m2 + (mVF + ^/^^WTW 

^3 = ^ . 1^ , (2.9) 

[2M2F + 8m2 + 2MVM2F2 + Am?F] ' 

with their propagators given by 



Gi = ^, (2.10) 



respectively, where 

m = ^ rnf + ml . (2-13) 

The propagators of the massive active-like mode and unsterile-like mode follow closely the 
propagators of the active-like mode and unsterile-like mode of Ref. |32j . 
As in Ref. [32], we assume that mi, 2 ^ M and self-consistently 
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m 



M2F(P2,3 

In this approximation, it is convenient to introduce 



<1. (2.14) 



ei,2(p) = -/^ ; e{p) = Jej{p) + el{p), (2.15) 
M \/F{p) V 



and the mixing angles 



„ ei mi . „ €2 m2 
cos u = — = — ; sm 6' = — = — 
em em 



1 e(v) 
cos$ = ; sin $ = . (2.16) 



We note that whereas the angle 9 does not depend on momentum, the angle $ does depend 
on momentum through -F(p). 
To leading order in e's we find 



M2{p) = -Me\p) ■ M^{p) = M[l + e\p)\ 



(2.17) 



and the fields associated with the propagating eigenstates given by eqs. (2.7) - (2.9) can be 
written as 



i^i = — sin 6 + cos 6 v^^ , 
vi = cos $ cos 9 + cos $ siu z/qj — sin $ ujj , 
= sin $ cos 9 + sin $ sin 6 Ua,^ + cos $ vu ■ 



(2.18) 
(2.19) 
(2.20) 



The properties of the propagators can be found by following the calculation in Ref. [32] 
and in the following, we will summarize the results. For details, see 



A. Poles and Spectral Density for the Massive Active-like Mode 



Near the pole, the propagator behaves as 

1 



2 ' 



where 



and 



with 



Mi 

M2 



Ml 



A 



MA 



m 
M2 



M2 

A2" 



(2.21) 



(2.22) 



(2.23) 



(2.24) 



This is an isolated pole below the unparticle threshold at = M^. This pole lies on the 
real axis and describes a massive stable active-like propagating mode. 

The massive active-like propagator also features an inherited spectral density 



P2[X) 



e(x) ^ x^" sin(27rr7) 



TT 



x + 1- ^x'^'^ cos(27r?7) + ^ x^'? sin(27r?7) 



(2.25) 
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where 



X 



M2 



(2.26) 



The non-vanishing spectral density above the unparticle threshold at p"^ = will lead 
to interesting new phenomena in the disappearance and appearance probabilities discussed 
below. 

B. Complex Poles and Spectral Density for the Unsterile-like Mode 



The solution describes a pole in the complex plane (a resonance) and it exists only for 



Near this pole we find 



with 



Re(a;) > , < r] < 1/3. 



p'^-Mi{p) p-^-Ml + iM^T 



1 — f] 



(2.27) 
(2.28) 



and 



Ml = 



r 



_}_ ( TXT] 
1 + Ai-'' COS 

V 1 — ?7 
7r?7 



M2 1 

-Ai-'' sin 



(2.29) 



(2.30) 



A^3 ^'^W-ri, 

The imaginary part is a consequence of the fact that the real part of the pole is above the 
unparticle continuum and it describes the decay of the unsterile like mode into the massive 
active-like mode and particles in the "hidden" conformal sector. This can be understood 
using a renormalization group resummation argument presented in Ref. [32]. The spectral 
density is obtained from the discontinuity across the real axis in the complex p^ plane 

, , Bfx) A x'' sin(7r?7) 



TT 



+ 



A x"^ sin(7rr7) 



2 • 



(2.31) 



X — A x'' cos(7rr7) 

For consistency of the see-saw mechanism, to give small masses to the active-like neutrinos 
it is required that 



A 



thus ensuring that M2 = ^MA < M. 



m 
M2 



M2 
A2" 



< 1, 



(2.32) 



III. TOWARDS UNDERSTANDING ACTIVE- ACTIVE OSCILLATION DYNAM- 
ICS. 



The familiar quantum mechanical description of neutrino oscillation invokes single par- 
ticle quantum mechanical states. The essential ingredients are the expectation values 



-iHt 



Vi). In quantum field theory, the (fermionic) fields are expanded in terms of 



creation and annihilation operators of quanta associated with a Fock basis, which is deter- 
mined by the spinor basis functions. For the mass eigenstates these are the solutions of 
the (free) Dirac equation, whereas for fiavor states such choice is ambiguous To focus 
the discussion on a simple case to be generalized later, consider a fiavor doublet of Dirac 
neutrinos with an off-diagonal mass term, with Lagrangian density 



£ = ^ 



M 



where 



(3.1) 



(3.2) 



Diagonalizing the mass matrix via a unitary transformation leads to the fields ipi^2 associated 
with the mass eigenstates mi^2, related to the fiavor fields ipe,n as 



ipe = COS Oipi + sin 9 ip2 , 
ip^i = cos 94'2 — sin 9 ipi ■ 



(3.3) 
(3.4) 



The neutrino fields '?/'i_2 are quantized as usual^ 

1 



^i(x,0) 



k,h=±l 



^|Jiik,0) = \h{k,h)Ui{k,h) + dl{k,h)Vi{-k,h) 



1,2 



(3.5) 



where the annihilation operators 6j, di and creation operators foj, obey the usual anticom- 
mutation relations and the spinors U, V are solutions of the Dirac equations 



[a-k + Pm4Ui{k,h) = Ei{k)U,{k,h), 
[a-k + l3mi\Vi{-k,h) = -Ei{k)Vi{-k,h) 



(3.6) 
(3.7) 



^ We keep the same notation for the field and its spatial Fourier transform to simplify notation, no confusion 
should arise since the arguments are different. 
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with Ei{k) = a/A;^ + mf, and are eigenstates of the hehcity operator with eigenvalues h = 
±1. Obviously there is no unambiguous quantization of the flavor spinor fields because 
there is no unambiguous choice of basis spinors U, V for these |35] . One could chose spinor 
solutions of the massless Dirac equation or with the masses in the diagonal entries in the 
mass matrix M. Any of these choices correspond to spinor solutions which are not orthogonal 
to any of the solutions of the mass eigenstates (either positive or negative energy) 
The positive energy single particle quantum mass eigenstates of helicity h are 

|z/,(fc, h)) = h\{k, h) |0) = 4{k, 0) |0) Ui{k, h) , (3.8) 

{iy,{k,h)\ = {0\hik,h) = Ulik,h){0\ilji{k,0). (3.9) 

We define the single particle flavor states as 

l^e) = cos6' |z/i) + sin^^ |z/2) , (3.10) 
lUfj) = cos^^ |z/2) — sin6' |z/i) , (3.11) 

where the quantum numbers k, h are common to all the states. Since the Hamiltonian is 
diagonal in the mass basis, the transition amplitudes (z/q| e~*^* (a,/3 = e,/i) require the 
overlaps 

{u,{k, /i)| e-^^* \u,{k, h)) = Ujik, h){0\ ^,{k, t)^,{k, 0) |0) 7°f/,(^, h) , (3.12) 

where ilJi{k,t) are the Heisenberg field operators. In terms of the propagator for t > 0, the 
overlaps are given by 

{u,{k, /i)| e-^^* \u,{k, h)) = Ujik, h) {iS^ik, t > 0))j'Ui{k, h) , (3.13) 

thereby establishing a direct relation between the quantum field theory propagators for the 
mass eigenstates and the single particle transition amplitudes. 
It is convenient to work in the chiral representation 

: '] ; .'-('°]. (3.14) 

in which we find the positive energy spinors for the mass eigenstates with mass rrii to be 



U,^,,n)^jm±±i I, (3,5) 



2k 



f/.A-l) = I . (3,16) 



2k 



v{k,-l) 



where Ei{k) = \/k'^ + mf and the two component spinors are eigenstates of hehcity 

a ■kv{k,h) = hv{k,h) ] h = ±1. (3.17) 

The spinors U are normahzed to unity and become right handed {h = 1) or left handed 
[h = —1) respectively in the ultrarelativistic limit. The propagator in real time is given by 

S^ik,t)= I ^e-^°*S.(p), (3.18) 



and for free field theory 



S^iP) = ^^^^^■ (3.19) 



It is straightforward to confirm that for t > the result of Eq. (3.13) is e~*-^»('^)*. 

In generalizing the above discussion to the case of active-unsterile mixing, we face two 
caveats: 



1. the mixing angle $ in Eq. (2.16) depends on the four momentum; this is one of the 
features that prompted the study of unsterile-active mixing as a potential reconciliation 
of MiniBoone and LSND data as suggested in Ref. [33j, 

2. the pole of the unsterile-like eigenstate z/3 is complex with the imaginary part being 
associated with the decay of the sterile neutrino into an active one and quanta of the 
hidden conformal sector 



A full quantum field theory study along the lines presented in Refs. [36l |37] incorporates 
the full propagators and spectral densities, which are dominated by the poles below the con- 
tinuum threshold in the case of the active-like modes and the complex pole above threshold 
for the unsterile like mode. The (four) momentum dependence of the mixing angles is then 
evaluated at the position of these poles (see Refs. [36], [37] for a detailed discussion). In these 
references, the dynamics of neutrino oscillations was studied as an initial value problem for 
an initial wave packet. However, the real oscillation experiment deals with neutrinos pro- 
duced at a source and detected via charged leptons in a far detector, so that the neutrino is 
virtual in this process. In this case, the real time quantum field theoretic analysis of oscil- 
lation dynamics is more complex. The interaction vertices are written in terms of the fields 
that create and annihilate the mass eigenstates and "flavor" neutrinos are an intermediate 
state, described by a propagator associated with the mass eigenstates. Such study is just 
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beginning to emerge, and progress will be reported elsewhere [38]. A full quantum field 
theoretical study of active-unsterile dynamics is beyond the scope of this article, but we can 
extrapolate from the study in Refs. [511 EI] take the mixing angles to be associated 
with the single (quasi) particle states, therefore evaluated on the mass shell of the mass 
eigenstates corresponding to the poles in the respective propagators. 

In the case of the unsterile mode, the imaginary part of the pole (decay width) is sublead- 



ing in the small ratio A (see Eq. (2.30)) and we neglect it in the evaluation of the mixing 
angle $ for the mass eigenstate z/3. 



Inverting the relation (2.18 2.20) and evaluating the mixing angles on the single particle 
mass shells we obtain, 



Wai) = COS $2 COS 6'|z/2) + sin $3 cos6'|i/3) — sin6'|z/i) , 
Wa2) = cos<l'2sin6'|z/2) + sin<l>3sin6'|z/3) + cos6'|t^i) , 
luu) = - sin $2 1 1^2) + cos $3 1 z/3) , 



(3.20) 
(3.21) 
(3.22) 



where to leading order in A we find 



sin $2 



A 



2[1 + 



-1 
2 J 



sin $,q 



5 



cos $5 



COS 



1 + 



and 



Ai-'' cos 



1 — 1] 



(3.23) 



(3.24) 



(3.25) 



We will focus on negative helicity positive energy states, for which the spinor wave functions 



are given by Eq. (3.16). 



An immediate caveat of this formulation is that as a consequence of the four-momentum 
dependence of the mixing angle $, the states i'ai,a2,u introduced above as a straightforward 
generalization of the familiar quantum mechanical description, are not orthogonal, despite 
the orthogonality of the mass eigenstates 1/1,2,3- This is because $2 7^ *&3- This is an 
unavoidable consequence of the energy-momentum dependence of the mixing angle and of 
the effort to establish a correspondence with the familiar single particle quantum mechanical 
description of an inherently many particle problem. The non-orthogonality of these states 
is small for A ^ 1 as manifest by Eqs. ( 3.23[3.24[3.25 ). 
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While the single particle quantum mechanical analogy has limitations in absence of a 
full quantum field theory treatment of neutrino mixing and oscillations directly in real time, 
we adopt the approximate single particle description afforded by ( 3.20[3.21[3.22 ) as a proxy 
description and proceed to explore this formulation as a prelude towards a firmer under- 
standing. A thorough field theoretical description of real time oscillations is postponed for 
future work. 

The fields associated with the mass eigenstates V'j, z = 1, 2, 3 are expanded at t = just 



as in Eq. (|3.5|) but with the spinors U, V solutions of the Dirac equations with masses 

TTrj ^ 



1 



Mi = 0, M2 = -MA, M3 = M 



1 + Ai-*) cos 



1 — r// . 



respectively |39]. The annihilation operators b{k, h), d{k, h) and creation operators 
b^{k,h), d'^{k,h) are interpolating Heisenberg field operators whose time evolution is de- 
termined by the total Hamiltonian such that [39] 



(0|^4A;,t)V'i(A;,0)|0) = ^i(A;,t>0). 



(3.26) 



In order to carry out the integration in po as in Eq. (3.18), it is convenient to write the 
propagators in terms of their spectral representation 



Jo 



p'^-Q^ + iO+ 



(3.27) 



from which we find 

Si{k,t > 0) 

where 



dQ' 



'i£{Q)t 



(3.28) 



(3.29) 



To leading order in A, we have pi^\Q^) = pf\Q'^) ; i = 2,3. For the massive active-like 
mode 



p^^\Q') = Z2 6{Q'-Ml) + 



P2{x) 
M2 



X 



M2 



(3.30) 



where Z2 is given by Eq. (2.22) and P2{x) by Eq. (2.25). For the unsterile like mode 



(3.31) 
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^ A=0.1; Ti=0.1; k=M 



0.8 



(Mt)^"^ I[Mt] 




Mt 



FIG. 1: The continuum integral I[Mt] = \{i^2\ e"*-^* |i^2>|contM2- 



where P3{x) is given by Eq. (2.31 ). Near the complex pole in the continuum, the propagator 



for the unsterile-like neutrino can be approximated the Breit-Wigner form Eq. (2.27), which 
for A,ri describes a narrow resonance near the real axis. 

For long times, the integral over the dispersive variable is dominated by the pole and the 
continuum threshold. The technical details of the evaluation of the threshold contribution 



are found in the appendix (see Eq. (A16)). For negative helicity states, we find 



(^1 












(3.32) 












\ ~iEM{k)t 


(3.33) 




^MHJ 


^2 




^2) 




rt 




(3.34) 



where Ei{k) = a/A;^ + Mf, T is given by Eq. (2.30), 7 = E2,{k)/M2 is the Lorentz time 
dilation factor and (see appendix, Eq. ( A17[ )) 

A2 



A. 



An 



sin(27rr/) r(l + 2r]) . 



(3.35) 



The extra factor 1/2 in (3.33) as compared to the bosonic case studied in the appendix (see 



Eq. (A16)) is a result of the spinor overlaps in the limit A <^ 1. 



The result in Eq. (3.33) is valid in the long time limit when the threshold contribution 
dominates the integral over the continuum |l0]. Fig. ([T]) displays the integral over the 
continuum and confirms its asymptotic long time limit. 
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We are now in position to obtain the active neutrino disappearance and appearance 
probabilities of active-like neutrinos. 



A. Disappearance probabilities 



The disappearance transition amplitudes are given by 



rai/ 



2 a,,. l^-i-f^< l^y^) +cOs2<l>2COs2e(z/2|e-*^*|z/2) 



sm U\Vi \e 

,2 „„„2 a /,, l-iHt 



+ sin^ $3 cos^ d (z/3 le"*""* , (3.36) 
cos^ 6{ui |e"*-^* \ui) + cos^ $2 sin^ 6 (1/3 \e~^^^ \v2) 
+ sin^ $3 sin^ 9 (1/3 je"^-^* lz/3) , (3.37) 



where the transition amplitudes \e |z/j) ;z = 1,2,3 are given by Eqs.(3.32 3.34). In 



obtaining the probabilities there are several interference terms which manifest themselves on 
different time scales. We will only keep the slow oscillatory terms involving the differences 

E^{k) - E2{k) ^ Ml/2k and Em(A;) - E^ik) ^ {Ml - M^)/2EM{k) 

in the ultrarelativistic limit, and neglect phases that oscillate on much shorter time scales, 
since these average out. As is customary, we replace time by the baseline t — )■ L, and define 

i+2-n 



A{L) 



Euik) 



(3.38) 



We then have 



Pai^ai (L) = Zl cos^ $2 cos^ e + sin^ ^ + — cos^ $2 sin^ IB cos 



MIL 



+ 2 A{L) Z3 cos^ $2 sin^ $3 cos* Q e cos 



+ Zl sin* $3 cos* e e 



2k 

{M^-Ml)L 
2EM{k) 



Va^^aiiL) = Z| cos^ $2 sin* ^ + COS* ^ + ^ cos^ $2 sin^ 2^ COS 



MIL 



+ 2 A{L) Z3 cos^ $2 sin^ $3 sin* 9e cos 
+ Z|sin*$3sin*^e~7^. 



2k 

{M^-Ml)L 
2EM{k) 



+ i2(L)cos*$2Cos*^ 

+ i 

(3.39) 

+ i2(2,)cos^<l>2sin*^ 

+ e 

(3.40) 
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B. Appearance probability 



Prom the transition amplitude 



|e-^^* \va2) = \ Sin cos^ $2 {^2 le"*^* ji^s) + sin^ $3 {u^ je"*^* I1/3) - {1^1 le"*^* } , 

(3.41) 

by keeping only the slow interference terms and again replacing t ^ L, where L is the 
baseline, we find 



sin^ 26 



1 + 



Zl + i^(L)J cos'' $2 - 2Z2 cos^ $2 cos 

(M2 - Ml)L 



MIL 
2k 



+ 2 A(L) Z3 cos^ $2 sin^ $36 "^f cos 



2£;m(A;) 



+ Z|sin^ ^ae'^-^ 



(3.42) 



where for A <S 1, it follows that 



Ml-M^^ M^A^ cos 



(3.43) 



C. Consequences of the anomalous dimension 



The unsterile neutrino is characterized by the anomalous dimension 77, which is responsible 
for A, Z2, Z^jT 7^ and $3 7^ $2- Therefore it is important to quantify the potential 
phenomenological effects of the appearance and disappearance probabilities associated with 
a non- vanishing (and perhaps large) anomalous dimension. For 77 = 0, these transition 
probabilites are the following (we neglect terms that oscillate on the time scale 1/M <^k/M.\ 
and average out) 



( cos^ $ cos^ 6 + sin^ 6*)^ + sin^ $ cos^ 9 - cos^ $ sin^(26') sin^ 



ai— )-ai 



MjL 

4k 



(3.44) 



( cos^ $ sin^ 6 + cos^ 6)^ + sin^ $ sin^ 6 - cos^ $ sin^ 26 sin^ 



MIL 
Ak J 



(3.45) 



ai— >-02 ' 02— )-<ii 



. 9 ^ ^ I sin'' $ 9 . 9 
sm^ 2^ < h cos $ sm^ 



MIL 
4A; J f ' 



(3.46) 
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where 

TTl 

sin $ = , . (3.47) 

In the above expressions, we have neglected interference terms that oscillate on the fast 
time scale oa 1/M which average out on the longer time scales of the oscillatory contributions 
displayed in these expressions^. 

We note that Vai^ai "^02-^^02 ■ This is a consequence of the fact that the transformation 
between flavor and mass eigenstates is not unitary in the active sector. In principle the 
existence of a canonical {rj = 0) sterile neutrino may be experimentally determined by 
measuring the disappearance probability for both active species (by measuring the associated 
charged leptons) along with the appearance probability. The difference in the disappearance 
probability for both active species signals the presence of a "sterile" degree of freedom that 
enters in the deflnition of the mass eigenstates but not in the weak interaction vertices. The 
combined measurement of all three probabilities for flxed baseline and energy would allow 
us to extract both mixing angles and oscillation lengths. 

A non-vanishing anomalous dimension introduces different angles $2,3 as a consequence 
of the momentum dependence of the mixing angle. This then gives rise to new oscillatory 
contribution with a different oscillation length that is also multiplied by an attenuation 
factor. To understand the different time scales it is convenient to introduce 

^ _ MIL _ M^A^L ^ ^^J^ 
2k Sk Lggc 

where Lgsc is the usual active oscillation length, as well as fl, the new oscillation frequency 
arising from the interference between the massive active-like threshold and the unsterile-like 
pole 



2k 



l-2t) 

A 1-" 



COS 



?7J 



in terms of which we flnd 

" " TTTj 



FL ~ 
— = \l tan 
27 



(3.49) 



(3.50) 



.I-77J 

For A <^ 1, the oscillation frequency Q, is larger than the active oscillation frequency Jl. As 
a consequence of this discrepancy in the frequency of the oscillatory contributions, as well 
as the attenuation and the wave- function renormalization factors, it follows that both the 



^ This is the reason that Va^^a2 does not vanish at t = 0. 
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FIG. 2: Disappearance and appearance probabilities for A = 0.1 and rj = 0,0.1 to compare 
canonical and unsterile neutrinos. 




FIG. 3: By the first oscillation peak, the difference between canonical sterile and unsterile neutrinos 
has disappeared. 



disappearance and appearance probabilities are suppressed relative to the case of a canonical 
sterile neutrino. However, because of the attenuation factors, the suppression is substantial 
only for short baseline experiments, namely ^ 1. This suppression is displayed in Fig. 
(|2|. We can also see in Fig. (|3| that by the first oscillation peak, the novel oscillatory 
contributions are well suppressed that it will be difficult to differentiate between canonical 
sterile and unsterile neutrinos. 

The discussion above applies in the limit A ^ 1 which has been invoked from the be- 
ginning to establish a see-saw hierarchy between the active-like and unsterile-like neutrinos. 
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FIG. 4: For larger A the suppression effect due to tlie anomalous dimension is more prominent. 



However, the proposal for a "solution" of the LSND/MiniBooNE discrepancy introduces a 
sterile neutrino in the eV mass range, namely within the same mass range as the active 
neutrinos. This possibility requires that A ~ Although our study does not apply 

directly to this regime, we can extrapolate to this regime in several aspects. 

First, the mixing angle $3, which determines the overlap between active-like and unsterile- 
like modes, becomes of the same order as the angle 6 that determines the overlap between 
active-like states. This modifies the appearance and disappearance probabilities by overall 
normalizations. However, the argument on Lorentz invariance still implies, at least within 
the quantum mechanical description of oscillations, that the angles are evaluated on the 
mass shell of the single particle states and do not depend on the energy. 



For A ~ 0{1), it follows from Eq. (3.49) that i7 ~ and the oscillation frequencies of 
both oscillatory terms in the probabilites are of the same order. Furthermore, the factor 



A of Eq. (3.38) that determines the interference terms from the threshold of the massive 
active mode and the unsterile mode become of order sin[27r?7]. Therefore, for rj < 1/3, 
but not too small, this overlap yields a potentially interesting energy dependence. In this 
case, the suppression of the probabilities as a consequence of the anomalous dimension is 
non- vanishing even for long baseline events, namely with Q ~ C^(l)- However, for large rj 



such that tcLrnvr] ~ 1 and A ~ 0{1), the attenuation length of the overlap (see Eq. (3.50)) 
is substantial and the interference between active and unsterile-like is suppressed over the 
baseline. The appearance and disappearance probabilities for A = 0.5 are shown in Fig. 
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IV. SUMMARY 



In this article, we considered the possibihty that the SU{2) singlet sterile neutrino is 
an unparticle and mixes with the active neutrinos. The unparticle sterile neutrino - or 
unsterile neutrino - is an interpolating field that describes a multiparticle continuum as 
a consequence of coupling to a "hidden" conformal sector and whose correlation functions 
feature an anomalous scaling dimension rj. We analyzed the consequences of its mixing with 
two active neutrinos via a "minimal" see-saw mass matrix, with massless active neutrinos 
and no active-active mixing term. 

We have introduced a generalization of the usual quantum mechanical description of the 
dynamics of mixing and oscillation that incorporates the propagators for the fields associated 
with the mass eigenstates and includes off-shell corrections to the time dependence of the ap- 
pearance and disappearance probabilities for the active-like neutrinos. We find a remarkable 
interference phenomenon between the massive active-like and unsterile-like modes as a con- 
sequence of threshold effects which modify both appearance and disappearance probabilities 
in a novel manner. 

The presence of the anomalous dimension that gives the unparticle nature to the ster- 
ile neutrino has profound consequences in the disappearance and appearance probabilities 
for the active-like modes. For a canonical sterile neutrino, the disappearance probability 
is different for the different fiavors as a consequence of the fact that the transformation 
between fiavor and mass eigenstates is not unitary in the active sector. This same feature 
remains in case of the unsterile neutrino, however, the anomalous dimension is responsible 
for novel time dependent phenomena corresponding to the new oscillatory term arising from 
the interference of the threshold of the massive active-like mode and the unsterile mode. 
This oscillatory term is multiplied by an attenuating function of baseline but the oscillation 
length is different from the active-active one. These novel contributions are consequences of 
the unparticle nature of the sterile neutrino and result in a suppression of both probabili- 
ties on short baseline experiments as compared to the case of active-sterile mixing with a 
canonical sterile neutrino. Combined analysis of short and long baseline experiments may, 
therefore, provide a diagnostic tool for the unparticle nature of a sterile neutrino. Another 
important manifestation of the unparticle nature of the sterile neutrino is the non-vanishing 
spectral density "inherited" by the massive active mode, which may lead to the possibility 
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of new reaction channels when the standard model interactions are present. 

Although these effects would be potentially relevant in the phenomenological reconcili- 
ation between the LSND and MiniBooNE data, we find that as a consequence of Lorentz 
invariance, the mixing angles, when evaluated on the mass shell of the corresponding mass 
eigenstates do not depend on energy as required for a reconciliation of the data[33]. Fur- 
thermore, within the simple model considered here, we find that either with a large see-saw 
hierarchy of masses between the unsterile-like and the active like neutrinos, or even if they 
all feature the same mass scale, the new interference terms that could be responsible for 
novel effects are probably too small to reconcile the LSND and MiniBooNE data. 

We would like to emphasize that there are, however, some important caveats: the frame- 
work introduced as a generalization of the quantum mechanical description of the dynamics 
of mixing and oscillations can at best be a proxy for a full quantum field theoretical treatment 
of the appearance and disappearance probabilities in real time in which "flavor" neutrinos 
are intermediate states, including the propagation of the unsterile degree of freedom. Fur- 
thermore, considering either a more realistic 3 + 1 scenario or a more general mixing matrix 
in 2 + 1 scheme with mass scales of the same order may provide a scenario in which the novel 
phenomena found here may yield rich phenomenology. For example, in 2 + 1 scheme, giving 
one of the active neutrinos a mass, which can be acquired by including another type beyond 
standard model sector, already results in a momentum dependent active-active mixing angle. 

Although the simple model studied here may not provide the reconciliation between the 
LSND and MiniBooNE data, the novel time dependent phenomena that emerges as a con- 
sequence of the unparticle nature of the sterile neutrino warrants further and deeper study. 
We anticipate important cosmological consequences in the equilibration of active neutrinos, 
novel mechanisms for sterile neutrino production and possibly interesting consequences for 
light dark matter candidates. 
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Appendix A: Oscillations of bosonic single particle states 



Neutrino mixing and oscillations is conventionally studied via a single particle quantum 
mechanical description. In this formulation, there is no distinction on whether the states 
of the quantum field theory correspond to fermionic or bosonic particles, and in the case of 
fermions, there is no mention of helicity or any other quantum numbers besides energy and 
momentum. 

In this appendix, by considering bosonic particles, we study the dynamics of oscillations to 
establish potential differences with the fermionic case described in the text. Let us consider 
a mixing of two canonical neutrinos. The mass eigenstates are given by 

Ui = cos^z/a + sin^z/f,, (Al) 
z/2 = cos — sin 9 Va , (A2) 

with z/a fe are the flavor eigenstates and the masses corresponding to the fields 1^12 are Mi 2, 
respectively. The transition amplitude |(t'a(^)|^'fe(0))P is calculated using ordinary quantum 
mechanics: 

Mt)\u,{0)) = {U,{t)\u,{t)) (Z^l(t)|z/l(0)) {M0)\u,{0)) + {Ua{t)\u2{t)) (z/2(t)k2(0)) (//2 (0) | /^^(O) ) , 

(A3) 



where the equal time overlaps (t'a,b|^i,2) are read from Eqs. (Al) and (A2), while the time 
evolution (2^1,2(^)1^1,2(0)) can be obtained from the field theory by considering the time 
evolution of a single particle state |lk)- To do that, let us first express the field describing 
the mass eigenstates using the creation and annihilation operators 

(k, t = 0) = (ak + aU) , (A4) 



where 



and E] = e + Mf. 



zy.(x, t = 0)= I ujik) e^'^-", (A5) 



(27r) 
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The time evolution (z/i,2(^) 1^^1,2(0)) is then given by 

{uM^M) = (lk|e-^*|lk) 

= 2Ej (0|z/j(k,t)z/](k,t = 0) |0) 



^E, I ^e--*G,(E,k) (A6) 



2E, 



dE e 



-iEt 



27rz E^-k^- Mf + ie 
e-'^^' . (A7) 



Therefore the amphtude is 

\{ua{t)\ui,{0))\'' = sin"2^^ sin 



2 _ „;„2ofl „;„2 iE2 - Ei) t 

2 



\{ML)\u,m\' ^ sin^2^ sin^ ^f'^^ , (A8) 

where on the second hne we have used the ultra-relativistic approximation. Here, L is the 
length of the baseline used in the experiment. 

Let us apply this to the 2 + 1-scenario. First of all, the overlaps (z/^^ 2 1^^1,2,3) can be derived 



from Eqs. (2.18) - (2.20). Let us start by expressing the energy eigenstates in terms of the 



"flavor" eigenstates as 

Wl)E,k = -Sm9E,kWai)k + COs9E,kWa2)k, (A9) 
W2)E,k = COS 4>E,k COS 9E,k Wai) k + COS 4>E,k Sm9E,k Wa2)k " siu Wu) k , (AlO) 

W3)E,k = Sm(l)E,k COs9E,k Wai)k + Sm(j)E,k Sm9E,k Wa2)k + COS(t)E,k \vu)k , (All) 

where we have made explicit the energy and momentum dependence. Since we are using 
one-particle quantum mechanics to describe the system, the energy eigenstates must be 
evaluated on-shell. Therefore 

\ui) = — sin I z/flj ) + cos 6* I z/a2 ) , (A12) 

W2) = cos$2 cos6' |z/ai) + cos$2 sin6' Iz/qj) — sin$2 l^^c/) , (A13) 
1^3) = sin$3 cos6' |z/aj) + sin$3 sin6' Iz/qj) + cos$3 |z/f/) , (A14) 

where 0j means that the angle (j) is evaluated at = Ai], and we have dropped the index 
for the angle 9 as it is a constant and does not depend on the four-momentum. We note 
that sin 03 and cos 03 are complex as A^3 > M. Here, we have also dropped the indices on 
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the states to simplify the notation. We can then solve Eqs. (A12) - (A14) for the "flavor" 
eigenstates to obtain the overlaps. 

Next, we need the time evolution of the mass eigenstates (i/j(t)|z/j(0)), especially for 
i = 2,3. We can obtain these by substituting their respective propagators into Eq. (A6) 



and it is convenient to use the dispersive form of the propagators [39]. The time evolution 
is then given by 



Z2 e-'^^^ + e"*-^^^* / dx p2{x) exp 



E 



M 



(A15) 



where we have introduced Euik) = \/k^ + M^, and x is the dimensionless variable defined 



in Eq. (2.24). 



For the case of z/2, since the pole is below the continuum, we can estimate the large time 
behavior of 5*2 by replacing p2 with its near-threshold behavior. Therefore 



A 2 roo 

Z2 e-'-^2(^)* + — sin(27rr7) e-'^^"^''^' / dx x^'' 
4vr Jo 



exp 



M2 ■ 
-i -^—tx 



E 



M 



Z2 e-'^'' + 



E 



M 



-iEMt 



with 



A2 = — sin(27r77) r(l + 2r/) . 
47r 



(A16) 



(A17) 



For z/3, we can estimate 6*3 by approximating the spectral density by a Breit-Wigner 



Lorentzian, (see Eq. (2.27)) from which we obtain 



27 



(A18) 



where F is the decay width given by Eq. (2.30) and 7 = E^{k)/M.-i is the Lorentz factor. 



23 



Finally, by combining all the results above, we find that at large time 



IK(t)|^a,(0))P 



sin^ 26 { 



^ |l + COS^<l>2 



^ X 2+4r, 



Zi + Ai 



+ 2 ^2 ^3 cos^ $2 I sin$3|^e 27* 

+ Z||sin$3|^e 
sin^ 26 



t ( Em 



— 2 Zo cos $2 cos 

cos [{Es - EM)t - "^{1 + 2r])] 

(A19) 



1 + cos $2 



2+iri 



— 2 Z2 cos^ $2 cos 



2k 



2 -hL 



+ 2yl2 ^3 cos $2 I sin$3| e 27 



E 



M 



COS 



{Ml-M')L TT 



2k 



(1 + 2//) 



Z| I sin $3|^e c , 



(A20) 



where in going from Eq. (A19) to Eq. (A20) we have used the ultra-relativistic approxima- 



tion and replaced t by the baseline L. 
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